Theory of helical spin crystals: phases, textures and properties 
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Motivated by recent experiments on the itinerant helimagnet MnSi, we study general properties of 
helical spin crystals - magnetic structures obtained by superposing distinct spin spirals. An effective 
Landau description of helical spin crystals is introduced and simple rules for stabilizing various spin 
crystal structures over single spirals are established. Curious properties of the magnetic structures 
so obtained, such as symmetry stabilized topological textures and missing Bragg reflections are 
pointed out. The response of helical spin crystals to crystalline anisotropy, magnetic field and non- 
magnetic disorder are studied, with special reference to the bccl spin structure, a promising starting 
point for discussing the 'partial order' phases seen at high pressure in MnSi. Similar approaches 
may be applied to other crystallization problems such as Larkin-Ovchinnikov-Fulde-Ferrel states in 
spin-imbalanced superconductors. 
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I. INTRODUCTION 

Long period helical spin order resulting from 
Dzyaloshinskii-Moriya (DM) spin-orbit coupling^ in 
non-centrosymmetric itinerant magnets (e.g. MnSi, 
FezCoi-zSi, FeGe with crystal structure B20) has been 
intensively studied in the pasti2iMiSiSiZi2i2ii2i One mate- 
rial of this family, MnSi, has attracted recent interest 
because of its puzzling behavior under applied hydro- 
static pressure! 11 ' 12 A state with non-Fermi liquid trans- 
port properties 1 ^ is obtained over a wide range of pres- 
sures above a critical threshold p c . Beginning at the 
same critical pressure, but over a smaller pressure range, 
magnetic 'partial order' is observed in neutron scattering 
experiments! 14 ' 15 While usual helical order gives rise to 
sharp Bragg peaks in neutron scattering (corresponding 
to the periodicity of the helical spin-density wave), 'par- 
tial order' is characterized by a neutron scattering signal 
which is smeared over a wavevector sphere rather than 
localized at discrete points in reciprocal space. 

Recent theoretical work on electronic properties, crit- 
ical fluctuations and collective modes of helimagnets 
are in Refs. 16 ' 17 ' 18 . Theoretical proposals for the high 
pressure state of MnSi have invoked proximity to a 
quantum multi-critical poinfeiS or magnetic liquid-gas 
transitions. 2 S Closest in spirit to our approach are the 
skyrmion-like magnetic patterns studied in Refsj 21 ' 22 . 

Recently, we have proposed a novel kind of magnetic 
order, the helical spin crystal, as a promising starting 
point for a theory of 'partial order'. 23 Helical spin crystals 
are magnetic patterns, which are obtained by superposi- 
tion of several helical spin-density waves which propagate 
in different directions. There is a substabtial resemblance 
to multi-fc magnetic structures (also known as multiple-g 
or multiple spin density wave states) i 24 ' 25 ' 26 ' 27 ' 28 But in 
contrast to most other magnetic multi-fc systems, in he- 
lical spin crystals the ordering wavevectors are selected 
from an infinite number of degenerate modes lying on a 
sphere in reciprocal space - a process analogous to the 
crystallization of liquids. 

In this work, we present a detailed theory of such struc- 



tures. The stability, structure and distinctive properties 
of such states are described, and the consequences of cou- 
pling to non-magnetic disorder is discussed. 



The paper is organized as follows. First, we review 
the standard theory of hclimagnetism in Section [Til nn ~ 
ishing with a short remark about more general helical 
magnetic states. Then, the theory of helical spin crystals 
is developed. The requirements to energetically stabi- 
lize helical spin crystal states are investigated in Section 
Mil The analysis works in two directions. First, we es- 
tablish a phase diagram in terms of natural parameters 
which tune the interaction between helical modes and 
second, we give simple rules to construct model interac- 
tions which stabilize a large class of helical spin crystals. 
The remaining parts of the paper are dedicated to ex- 
tracting testable consequences of these novel magnetic 
states. In Section IIV1 we give a description of the most 
prominent spin crystals which emerge from our energetic 
analysis in terms of their symmetry. It is shown that 
the symmetry of the magnetic state may stabilize topo- 
logical textures like merons and anti-vortices which are 
otherwise not expected to be stable in the present con- 
text, given the order parameter and dimensionality of 
the system. Symmetry also determines, which higher- 
harmonics Bragg peaks these structures would produce. 
We subsequently study the response of helical spin crys- 
tals with respect to different perturbations in Section 
fVl For example, sub-leading spin orbit coupling (crys- 
tal anisotropy) locks the magnetic crystal to the under- 
lying atomic lattice and thus determines the location of 
magnetic Bragg peaks. We also study the response to 
an external magnetic field which, apart from producing 
a uniform magnetic moment, also leads to distinctive dis- 
tortions of the helical magnetic structure, which could be 
observable by neutron scattering. Finally, in Section IVII 
we investigate the implications of non-magnetic impuri- 
ties, which are expected to destroy long-range magnetic 
order and produce diffuse scattering. 
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II. LANDAU-GINZBURG THEORY OF 
HELIMAGNETISM 

For a cubic magnet without a center of inversion, the 
Landau-Ginzburg free energy to quadratic order in the 
magnetization M(r) is 

F 2 = (r M 2 + J(d a M p ){d a Mp) + 2DM ■ (V x M)> , 

(1) 

where (. . .) indicates sample averaging, r$,J,D are pa- 
rameters ( J > 0) and Einstein summation is understood. 
The last term of Eq. ([T]) is the DM interaction, which is 
odd under spatial inversion and originates in spin-orbit 
coupling^ Fourier transformation, M(r) = X) q m q eiq r 
with m_ q = m* , leads to 

F 2 =Y1 i( r + Jl 2 ) l™q| 2 + 2 ^ m q ' (Ml >< m q)] ■ ( 2 ) 

q 

Clearly, the energy is minimal for circularly polarized 
spiral modes, where V x M points in the direction of 
— 13 M. For such modes, 

F 2 = 5>( g )|m q | 2 , (3) 
q 

where r(q) = r - JQ 2 + J(q - Q) 2 with Q = \D\/J. 
The Gaussian theory thus determines both the chiral- 
ity of low-energy helical modes and their wavelength 
A = 2-n/Q. The latter is typically long (between 180A 
in MnSi and 2300A in Feo.3Coo.7Si), reflecting the small- 
ness of spin-orbit coupling effects compared to exchange. 
However, no preferred spiraling direction is selected by 
Eq. ([T]), since F2 is rotation- invariant. Cubic anisotropy 
terms which break this invariance are of higher order in 
the spin-orbit interaction and therefore small. We neglect 
them for the moment and reintroduce them later. 

The isotropic Gaussian theory leaves us with an infi- 
nite number of modes which become soft as r(Q) — > 0. 
They consist of helical spin-density waves with given chi- 
rality (determined by the sign of £>), whose wave- vectors 
lie on a sphere |q| = Q in reciprocal space. Each of 
these helical modes is determined by an amplitude and a 
phase. Hence, for each point q on the sphere, we define 
a complex order parameter i/> q (with ip—q = tpq) through 

mq= 2^l( ? q+ i? q)> ( 4 ) 

where e q , e q , and q, are mutually orthogonal unit vectors 
(with a defined handedness, given by the sign of D). Ob- 
viously, changing the phase of tpq is equivalent to rotating 
e q and e q around q. The phase of tpq is thus only defined 
relative to some initial choice of e q . The neutron scat- 
tering intensity is proportional to |q x m q | 2 = l/2|?/> q | 2 , 
independent of the phase. Changing the phase of tpq is 
also equivalent to translating M(r) along q. 

In the following, we study minima of the free energy 
in the ordered phase [r(Q) < 0]. These depend on the 



interactions between degenerate modes (i.e., free energy 
contributions, which are quartic or higher order in M). 
We only consider interactions which, as F2, have full ro- 
tation symmetry and we will include the weak crystal 
anisotropy last. The most general quartic term which 
has full rotation symmetry (transforming space and spin 
together) is of the form 

F *= C/ ( qi ' q2 ' ' m ^ ( m i 3 ' ' 

qi>q2,q3 
with q 4 = -(qi + q 2 + q 3 ). 

A. Single-spiral state 

For example, if f/(qi, q2, q3) is a constant, then F4 cx 
(M 4 ). If the interaction depends only on the local mag- 
netization amplitude, i.e., in general if F = i 7 2 + (/(M 2 )) 
for some function /, then the absolute minimum of F is 
given by a single-spiral state (also known as helical spin 
density wave) M(r) = mke lk ' r + m£e~ lk ' r , where a single 
pair of opposite momenta ±k is selected. To proof this, 
we write F as 

5>(g) - r(Q)] |m q | 2 + (r(Q)M 2 + /(M 2 )) . (6) 
q 

In the single-spiral state, M 2 is constant in space and 
it minimizes the first and the second term of Eq. (JB|) 
independently. Therefore, no other magnetic state can 
be lower in energy. 

Because Q is small, the relevant wavevectors entering 
Eq. ([5]) are also small and t/(qi, q 2 , qa) is effectively close 
to a constant. Therefore, the single-spiral state, as ob- 
served in Fe^Coi^Si, FeGe and in MnSi at ambient pres- 
sure, is the most natural helical magnetic order from the 
point of view of Landau theory. 

B. Linear superpositions of single-spiral states 

Motivated by the phenomenology of 'partial order', we 
will now extend the theory beyond this standard solution. 
We speculate that £/(qi, q 2 , q 3 ) is not constant, such that 
i*4 favors a linear superpositions of multiple spin-spirals 
with different wave-vectors on the sphere of degenerate 
modes |q| = Q. 

One may first speculate about magnetic patterns 
whose Fourier transform is non-zero everywhere on the 
wave- vector sphere and peaked infinitely sharply perpen- 
dicular to the sphere, i.e. 

I^ql 2 oc tf(|q| - Q) (7) 

[see Eq. Q]. This idea turns out to be complicated for 
at least two reasons. 

The first complication is that there is no continuous 
way of attributing a finite-amplitude spiral mode to each 
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point on the wavevector sphere. This can be seen by 
noting that e q [Eq. Q] is a tangent vector field on the 
sphere. Thus, it cannot be continuous (impossibility of 
combing a hedgehog)^ Thus there is no "uniform" su- 
perposition of helical modes on the sphere. The problem 
of singularities can be avoided if one assumes a ipq with 
point nodes on the sphere. 

The second complication is that higher harmonics 
would result in a broadening of the delta-function in 
Eq. 0. This is seen as follows. Consider three momenta 
qi , q2 , q3 on the wavevector sphere and q4 which is off 
the sphere. The non- vanishing modes m qi , m q2 , m q3 
couple linearly to m q4 via Eq. ([5]) and thus induce a 
higher harmonic "off-shell" mode m q4 ^ 0. Since this 
happens for every point away from the sphere, the ef- 
fect is an intrinsic broadening of the peak in |V> q | 2 , hi 
contradiction with the initial assumption of Eq. ([7]) . 



III. ENERGETICS OF HELICAL SPIN 
CRYSTALS 

In the following, we study magnetic structures which 
are superpositions of a finite number of degenerate heli- 
cal modes ipj with wavevectors ±kj, j = 1, . . . , N. We 
call the resulting states helical spin crystals, because of 
the analogy with weak crystallization theory of the solid- 
liquid transition^ 



A. Structure of the quartic interaction 

We assume that F4 is small, and that its main effect is 
to provide an interaction between the modes which are 
degenerate under F%. Thus, the relevant terms of F4 are 
those with |qi| = | q2 1 = |q3| = |q4 1 = Q- This phase- 
space constraint and rotational symmetry implies that 
the coupling function U depends only on two relative 
angles between the momenta 

£/( qi ,q 2 ,q 3 )|| qi | =Q = E/(M), (8) 
where we have chosen the following parameterization: 

29 = arccos(qi • q 2 ) 

r(q2-qi)-q3l ( 9 ) 

0/2 = arccos . 

1 - qi • q 2 

Geometrically, 0/2 is the angle between the two planes 
spanned by (qi,q 2 ) and (q3,q4) (Fig. [I]). In the special 
case qi + q 2 = 0, it becomes the angle between q 2 and 
q3. This mapping allows 9 and (f> to be interpreted as the 
polar and azimuthal angles of a sphere and the coupling 
U(9, </>) is a function on that sphere. Since it describes 
an effective coupling between modes on the wavevector 
sphere, the coupling U (9, 4>) has a status similar to that 
of Fermi liquid parameters in the theory of metals. 




FIG. 1: The set of quartets (qi, . - . , q4) satisfying |q;| = Q 
and qi + . . . + q4 = 0, modulo global rotations, may be pa- 
rameterized by two angles 9 and <j> as shown in this figure. 



The Landau free energy F = F2 + F4 of helical spin 
crystal states is calculated as follows. Obviously, 

F2 = r{Q)Y J \n 2 - (10) 

3 

The quartic term, Eq. © , may be split into three distinct 
contributions F4 = F s + F p + F nt . The first term, F s , is 
the self-interaction of each spiral mode with itself: 

F S = L/ S £>/, (11) 
3 

where U s = U(9 = 7r/2, cf> = 0). A minimum requirement 
for stability of the theory is U s > 0. 

The next term, F p , stems from pairwise interactions 
between modes. It is of the form 

F p = 2^y p (%)|^| 2 |^f, (12) 

i<3 

where 29 ij = arccos(kj • kj) and 

V p (9) ^U(^,4d)+sm 4 9U(9,0) + cos 4 9U(^~9,0). (13) 

Since \img^oV p (9) = 2U S is large and positive, F p pro- 
vides an efficient repulsion of modes which are too close 
on the wavevector sphere. Energetically stable super- 
positions of spirals require substantially smaller values 
of V p {9) and therefore big-enough angles between them 
( Section IIII Cp . This "mode repulsion" suggests that a 
continuous distribution of Fourier modes on the wavevec- 
tor sphere, as discussed in the Section III Bl is difficult to 
stabilize energetically by slowly varying coupling func- 
tions U(9, (j)). However, singular behavior of the coupling 
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function is conceivable when discussing effective interac- 
tions of low energy modest 

Finally, the non-trivial quartic term F nt stems from 
quartets of modes, whose wavevectors sum up to zero. 
The geometry of four equal-length wavevectors summing 
to zero is depicted in Fig. [TJ A quartet-contribution oc- 
curs for < 29 < tv and < cj)/2 < it. In this case, 
the eight wavevectors ±kj form the vertices of a cuboid. 
Hence, 

J ' — ' U 

where the summation is over such quartets. For ex- 
ample, if kfj + kj 2 + kj 3 + k J4 = 0, there is a term 
Fj 1 .j 2 .j 3 j 4 oc '4>j 1 '4>j 2 'ipj 3 4'jA- The algebra of these terms, 
which depend on the relative phases of modes, is rather 
lengthy and obviously depending on the phase convention 
used to define the ip- variables. 

The term Fat fixes the relative phases of ipj to minimize 
the energy. If there is no frustration between minimizing 
each individual quartet-term, the phases arrange in such 
a way that all Fj 1 j 2 j 3 j 4 < 0. It follows that F nt < 0, 
after minimization for the V'-phases. A special case is 
the quartet of modes, whose wavevectors ±kj form the 
vertices of a cube (e.g. modes along (111)). In this case, 
Fj 1 .j 2 .j 3 ,j i — independently of the coupling function 
U (9, </>) as a consequence of rotational symmetries^ This 
is the reason why no such quartet term appears in the 
theory of Bak and Jensen^ As a consequence, a superpo- 
sition of four modes along all (111) directions is necessar- 
ily unstable towards shifting the wavevectors away from 
the perfect cubic configuration in order to gain energy 
from F nt (see description of fee* below). 

B. Phase diagram 

If U (9, 4>) is only slowly varying, it is justified to expand 
it in spherical harmonics (Yi m ). That is, 

U(0, <p) = U + U n sin 9 cos 4> 

+ C7 20 (3 cos 2 9-1) + U 22 sin 2 9 cos 2(j), ^ 

where we retained all terms with I < 2, which satisfy the 
relation {7(0, 0) = U{n - 0, <£) = £7(0, 2tt - <j>). 

Different quartic expressions in terms of the real-space 
magnetization M(r) may lead to the same projected cou- 
pling U(9,(f>). For example, the three terms 

(Vm 4 + W [V (M 2 )] 2 + W" (d a )(d a Mp ) M 2 ^ 

(16) 

lead to E/(qi,q 2 ,q 3 ) = W + W'(qi +q 2 ) 2 - W"qi • q 2 - 
When projected onto the wavevector sphere, they gener- 
ate only two spherical harmonics, namely 

U = W + \{W + W")Q 2 : 

2 (17) 
U m = g(2W" - W")Q 2 . 
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FIG. 2: Phase diagram for W > and W' = 0. In the 
grey region, FA < and the quartic theory is unstable. The 
various phases are explained in the text. In contrast to our 
earlier use of these symbols,»23 "A" and "□" denote general 
states with 3 and 4 helical modes, respectively. 

The terms W and W have intuitive interpretations: W 
restricts the magnetization amplitude, whereas W fa- 
vors or disfavors modulations of M 2 , depending on the 
sign4£ In the following, we set W" = and use the four 
parameters W,W',Un and (7 22 to tune the interaction 
U{6,f). 

It is hard to find the exact global minimum of _F 2 + F4 
for a general coupling function 1/(9,4)). We therefore 
restrict ourselves to a certain variational class of mag- 
netic states and determine the minimum within this class. 
Recently^ we studied only states which can be obtained 
by superposition of those six modes which propagate 
along the (110) directions (6-mode model). Here, we ex- 
plore a much broader class of states. We include 

1. any superposition of those 13 spin-spirals which 
propagate along the directions 4 ^ (HI), (HO) and 
(100) (13 complex variables as variational pa- 
rameters) , 

2. any superposition of up to four spirals with arbi- 
trary propagation direction (four variables ipj and 
five independent angles between wavevectors as 
variational parameters) . 

Within these constraints, we have computed 45 a phase 
diagram as a function of the coupling parameters 
W, W, Uu and U22, shown in Figs. [3]and[! 

All magnetic ground states are equal-amplitude super- 
positions of 1,2,3,4 or 6 spiral modes. The body-centered 
cubic (bec) states are superpositions of all (llO)-modes. 
beel and bcc2 differ by the relative phases of the six 
interfering helical modes (see Section IIV A ip . The sim- 
ple cubic (sc) crystal consists of three mutually orthog- 
onal spirals (e.g. along all (100) directions). A face- 
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FIG. 3: Same as Fig. 1 with Q 2 W = -0.5W. The point 
(7n = U22 = is now at the phase boundary between spiral 
order and simple cubic. 
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FIG. 4: a) and b) For the locations A, . . . , E in the phase 
diagram of Fig. [2] the ratio between the pair interaction V p (0) 
and the self-interaction U a is plotted as a function of the angle 
29 (the angle between propagation directions of modes). A: 
single-spiral state; B: phase boundary between single-spiral 
and bccl; C: bccl; D: sc; E: A. 



centered cubic (fee) helical spin crystal is obtained by 
superposing all four (lll)-modes. However, the ground 
state is not fee but a small distortion of it: fee*. In fee*, 
the wavevectors are shifted slightly away from (111) in 
order to gain energy from the quartet-term F nt £* The 
symbols "A" and "□" are used differently here than in 



our previous paper^ Here, "□" stands for a superposi- 
tion of four modes with wavevectors as shown in Fig. [1] 
with 4>/2 = tt/2. Hence, wavevectors ±kj form a square 
cuboid and allow for one quartet-term F nt . The angle 29 
changes as a function of interaction parameters within 
the range 0.24tt < 29 < 0.38tt. Finally, the phase "A" 
consists of three modes. The wavevectors ki , k2, k3 point 
to the vertices of an equilateral triangle on the sphere, 
whose size is determined by the requirement that the 
mutual angle between two wavevectors, 29, minimizes 
Eq. (flU)) . The angle 29 is parameter-dependent and lies 
in the range 0.147T < 29 < 0.247T. 

As expected, a negative W [Eq. (fl"6)) ] favors multi- 
mode spin crystal states with varying magnetization am- 
plitude relative to the spiral state with constant M 2 
(compare Figs. and \3§ . Positive W has the opposite 
effect, and enhances the region of the spiral phase (not 
shown). The term W' alone (i.e., with Un = U22 = 0) 
stabilizes sc in the regime Q 2 W < —W/2 < 0. However, 
a small positive U22 is sufficient to favors bccl over sc. 

In conclusion, we observe that two helical spin crystals, 
bccl and sc, appear adjacent to the single-spiral state 
and are stable at relatively small values of Q 2 W', Un 
and 1/22- In the following, we study the properties of 
the bccl, and sc states since they are the most likely 
candidates of helical spin crystals from the point of view 
of energetics. 

C. Model interactions with exact ground states 

In the preceding Section, we established a variational 
phase diagram for "natural" , i.e., slowly varying coupling 
functions U(9,4>). Most phases in this phase diagram, 
("spiral", sc, bccl, bcc2, fee and "A") can be shown to 
be the exact global minima for some fine-tuned model 
interaction, that are constructed below. 

Let us consider the toy model F toy = F% + F s + F p 
[Eqs. (|10H12p ] where the quartet term F nt is dropped and 
we replace V p (9) by a constant V. In this model, local 
minima with N non- vanishing modes (N > 1) must have 
equal amplitudes \ipj\ 2 = \r\/(U s - V + NV)/2 and the 
minimum energy with N modes is 

1 r 2 

F t oy,N = -J ~ ~ ■ Vr _ v . (18) 
V + N 

There are three regimes. If < U s < V, the ground 
state is the single-mode phase with N = 1 (all but one ip 
are zero). For < V < U s , the energy is monotonically 
decreasing with N. It means that the system includes as 
many modes as possible to lower its energy. Finally, for 
V < or Us < 0, .Ftoy is not bounded from below and 
therefore unstable. 

To make use of this toy model, we must tune the in- 
teraction such that U(9, 4>) — > unless 9 = n/2 or <f) = 0. 
This removes the quartet term F nt and we are left with 
a model similar to F toy with the difference that V p (9) is 
not constant. 
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We now tune the interaction such that V p (0) is very 
kjg33 everywhere except at some angle 9 opt , where it 
has a narrow minimum. For such a pair interaction, 
all arrangements of modes which involve angles other 
than 2(9 opt are excluded. Within this constraint, we are 
left with our toy model with a constant pair interaction 
V = Vp(8 pt), but the number of modes is restricted to 
1 < N < 3, since no more than three modes can have 
equal mutual angles between them. It is clear that in the 
region — 3U s /2 < V p {9 opt ) < U s , the ground state is a 
helical spin crystal with three modes, i.e., the state "A" 
or sc (in the case 29 opt — 7r/2). 

A different class of exact ground states is obtained, if 
the interaction is tuned such that 

1. V p (9) is very big^ in the region < 29 < 29 c and 

2. V p (9) = V (constant) in the region 29 c < 29 < tt/2 

for some critical angle 29 c . In this way, modes whose 
wavevectors are too close are excluded. That is, it en- 
forces a "hard sphere" constraint |kj±kj| > 2(1 — cos 2 9 C ) 
on the wavevectors. The interaction between modes 
which satisfy this constraint is constant and reduces to 
the toy model. Therefore, if V < U s the ground state will 
include as many modes as geometrically possible by the 
"hard sphere" constraint. In the case 29 c = arccosl/3, 
we obtain fee, which has wavevectors at the vertices of 
a cube. In the case 29 c = 7r/3, we obtain bec whose 
wavevectors are the vertices of a cuboctahedron. This 
can be seen as follows. A real-space bec lattice corre- 
sponds to a fee reciprocal lattice and the wavevectors 
±kj of the bec spin crystals are the 12 nearest neighbors 
of the origin in the fee reciprocal lattice. Because fee 
is the cubic close packing of spheres, this represents the 
only arrangement of twelve vectors | ky | = Q which sat- 
isfies the constraint |kj — k,! > 1. (The hexagonal close 
packing is not acceptable because it does not consist of 
pairs ±kj.) In principle, this construction can be used to 
create models whose ground state contains an arbitrarily 
high number of modes. 

We can now use the insight from these constructed 
models to understand certain features of the phase dia- 
gram of Figs.[2]and[3] In Fig. [4] we plot V p {9) as obtained 
from the expansion Eq. (fT5j) at different places in the 
phase diagram of Fig. [5J In the beel region of the phase 
diagram (curve C), V p {9) is similar to the one constructed 
above: close to constant and small for 7r/3 < 29 < tt/2, 
big for 29 < tt/3. At the phase boundary between single- 
spiral and beel (curve B), V p is still bigger than U s , indi- 
cating that the quartet term F nt is essential to stabilize 
beel close to the phase boundary. The phase boundaries 
between single-spiral and three- mode states (sc or A) 
are exactly determined by the crossing of the minimum 
of V p (9) with U s , as illustrated by curves D and E. 



IV. STRUCTURE OF HELICAL SPIN 
CRYSTALS 

A. Symmetry properties and topology of helical 
spin crystals. 

Time-reversal symmetry (T) reverses the magnetiza- 
tion direction M — * — M and may be implemented in 
terms of the variables as ipj — > — ipj. 

Under spatial translation by a vector a, the ^-variables 
transform as ^ = ^ exp(ikja), i.e., they experience a 
phase change. 

If the set of vectors k, is linearly independent, every 
phase-change of modes simply amounts to a global trans- 
lation. This is the case in helical spin crystals with up to 
three non-planar modes, e.g. in sc and A. These states 
are periodic with Bravais vectors ai, a.2, a.3, which are de- 
termined by a.i ■ kj = 2nSij . It follows that time-reversal 
T is equivalent to a translation by (ai + &2 + 0-3) /2 for 
these states, which leads to a 2-fold symmetry inside the 
unit cell reminiscent of antiferromagnetism. 

In contrast, helical spin crystals with more than three 
modes like bec, fee and □ depend essentially on the rela- 
tive phases between the helical modes. For these states, 
T is not equivalent to any translation, since it is not 
possible that exp(ik,a) = —1 for all k 3 . As a conse- 
quence, these states are doubly degenerate in addition to 
the translational degeneracy. 

To obtain an understanding of the real-space picture 
of helical spin crystals, it is useful to study their sym- 
metries under rotations (reflections change the chirality 
and therefore never appear in the symmetry group). For 
example, if the structure M(r) has a n— fold rotation axis 
with direction u, then M || u along this axis. Therefore, 
Mi, which is the projection to the plane orthogonal to 
u, has a node at the axis. The winding number 2 - (vortic- 
ity) of the node is restricted by symmetry to the values 
1, 1 ± n, 1 ± 2n, etc. The resulting pattern in the vicin- 
ity of the rotation axis, M pointing along u in the center 
and Mi winding around it, resembles that of a skyrmion 
or meron. Such patterns are currently being discussed in 
the context of helimagnets i 2 1 ' 22 Here, we observe that 
they naturally appear resulting from rotational symme- 
tries. The simplest cases are winding numbers of +1 (for 
any n-fold axis) or —1 (only for n = 2). 

Apart from proper rotations, the point group may 
contain anti-rotations, i.e., symmetry operations which 
are composed of a rotation followed by T. n-fold anti- 
rotation axes are only possible for even numbers n, since 
they imply n/2-fold rotation symmetry. In the special 
case n — 2, it merely follows that M _L u along the axis, 
where u is the axis direction. Higher anti-rotation sym- 
metries with n = 4, 6, . . . imply M = along the axis, 
i.e., they create a line-node in the magnetization. In the 
vicinity of this line-node, M is approximately orthogonal 
to u (as can be seen by expanding M(r) to linear order 
around a point on the axis) . The winding number of Mi 
around the node, is restricted by symmetry to the values 
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FIG. 5: Geometry of the wavevectors ki,...,k6, which con- 
stitute the bcc spin crystal states. The vectors ±kj form the 
vertices of a cuboctahedron. 

1 ± n/2, 1 ± 3n/2, 1 ± 5n/2, etc. The simplest case for a 
4- fold anti- rotation axis is a winding number of —1 (anti- 
vortex line) and the simplest case for n = 6 is a winding 
number of —2. 

We have thus demonstrated the emergence of topolog- 
ical objects like merons and antivortices, which are not 
expected to be stable in the present context of a vectorial 
order parameter in three dimensions but which are stabi- 
lized by symmetry. Thus, rotation axes are meron lines 
and 4- fold anti-rotation axes are anti- vortex node-lines. 



1. Symmetry and real-space picture of bed 

The bcc helical spin crystals consist of six helical modes 
with wavevectors ki = [110], k 2 = [110], k 3 = [Oil], 
k 4 = [Oil], k 5 = [101] and k 6 = [101], as shown in 
Fig. and has the periodicity of a body centered cubic 
lattice. We chose the convention [see Eq. Q] 



for j = 1, . . . , 6 with z = [001], and we consider negative 
chirality [i.e., M • (V x M) < 0], such that e'j = k, x e". 
The geometry of the wavevectors allows for three quartet 
terms in the free energy produced by T x = ipl '^2'4 ! 5'4>&, 
T y = ^*^^3tp4, and T z — — V>3^|'05V'6- The sign in the 
definition of T z has been introduced for convenience. 

The transformation properties of the ^-variables and 
the three quartic terms under rotations are shown in Ta- 
ble |TJ From this, it can be deduced that the (rotation- 
invariant) quartet contribution to the free energy is 

F nt = 2A nt Re(T x +T y + T z ). (20) 

The value of the parameter A n t is not determined by sym- 
metry. Direct calculation yields 

A nt = ^[|,2arccos(i)]-if/(J^). (21) 

There are thus two bcc states, depending on the sign 
of A nt . In beel (bcc2), which corresponds to A nt > 
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i>3 


itp2 




T* 


T z 





TABLE I: Transformation properties of the ^-variables and 
three quartic terms (defined in Section TlV A 1|) of the bcc spin 
crystals under rotations. R z and R x , respectively, are 7r/2 
rotations around the z- and a;-axis. These two rotations gen- 
erate the cubic point group O and therefore, the behavior 
under any rotation which maps the 12 wavevectors onto each 
other may be obtained by combining these two operations. 



(< 0), the phases of the ^'s are such that T x , T y , T z are 
all negative (positive). Three out of the six phases are ar- 
bitrary due to global translation symmetry. This means 
that the magnetic pattern of beel and bcc2 is uniquely 
determined up to translational and time-reversal degen- 
eracy. 

The solution for A nt > 0, beel, turns out to be the bcc 
structure with the highest point group symmetry. By 
selecting the coordinate origin conveniently, we obtain 
— ipi = ip2 = —iip3 = iip4 = iips = —iip6 = SMq for beel, 
where S — ±1 is the time-reversal symmetry label and 
Mq > is the amplitude. From Table U we deduce that 
M(r) changes sign under a ir/2 rotation about the x, y 
or z axis. That is, the magnetic point group is 0(T) (in- 
ternational notation 432) with 4-fold anti-rotation axes 
at (100), 3-fold rotation axes at (111) and 2-fold anti- 
rotation axes at (110). 

The real-space representation of the beel state is 

(V2 s x (c y - c z ) - 2 s y s z \ 
V2s y {c z ~c x )-2s zSx , (22) 
V2 s z (c x - Cy) - 2 S x Sy j 

where s x — sin(Qx/v / 2), c x — cos(Qx/V2), etc. The re- 
sulting pattern was shown in Fig. 2 of our earlier paper^ 
In Fig. [6l we show the symmetry axes. As discussed 
above, the magnetization must vanish along the 4-fold 
anti-rotation axes, which are anti-vortices with winding 
number —1. The x, y, z axes, and their translations ac- 
cording to the bcc periodicity, form two interpenetrating 
cubic latices of such line-nodes. The cubic space diag- 
onals are 3-fold and the red arrowed lines of Fig. [S] are 
2-fold rotation axes. In the vicinity of these lines, the 
magnetization field is skyrmion-like (i.e., has a wind- 
ing number of +1). 

The fact that the beel state breaks T in a non-trivial 
way and cannot be restored by any translation is mani- 
fest in the occurrence of the T-breaking order parameter 
(M x M y M z ) = SM$/2 ^ 0, which is a magnetic octupole. 
This curious property may lead to distinctive anomalous 
effects, e.g. in the magnetotransport. 34 Octupolar mag- 
netic ordering has recently been discussed in different 
contexts^ 



8 



±5 



FIG. 6: (Color online). Symmetry of the bccl state. The fig- 
ure shows a cubic unit cell of bccl. Black lines are anti- vortex 
lines with 4-fold anti-rotation symmetry and vanishing mag- 
netization. The red (dark gray) lines are 2-fold rotation axes 
and the arrows indicate the direction of M. The structure 
has 3-fold rotation symmetry about all cubic space diagonals. 
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TABLE II: Transformation for the i/j-variables of sc under 
rotations. R z and R x are 7r/2-rotations, R\ux] is a 2%/ 3- 
rotation and iiriio] a 7r-rotation around the indicated axis. 



where s x — sin[Q(x + A/8)], c x — cos[Q(a- + A/8)], etc. 



B. Higher harmonics Fourier modes 



As briefly mentioned in Section lll Bl magnetic ordering 
in wavevectors ±kj generally induces higher harmonics 
in the magnetic structure. In the presence of magnetic 
order m^. ^ 0, the Landau free energy for the modes 
m q , which do not belong to the set , is (to quartic 
order) 



AF 



^>v,(q)|m q | 
q 



hv-(q) 



(24) 



with r^,(q) = r(q) + O (|0j| 2 ) and 
M<l) = -4 ^2 c/ (qi.q2,q3) (m qi • m q2 ) m q3 , (25) 



where the sum is restricted to qi,q2,q3 € {±kj} such 
that qi + q2 + q3 = q. The origin of the exchange field 
h^, is the coupling term Eq. |5]). In the following, we 
assume that f (q) > 0. Obviously, Eq. ([24]) then leads to 
induced modes 



m a = 



2. Symmetry and real-space picture of sc 

The simple cubic (sc) helical spin crystal consists of 
three modes with ki = [100], k 2 = [010] and k 3 = [001]. 
It forms a periodic structure with a cubic unit cell and 
the lattice constant is A = 2n/Q. The convention for e" 
is given by Eq. (fTTj|). where the unit vector z is replaced 
by [111]. 

The transformation properties of the tp- variables under 
rotations are given in Table [Hi By choosing the center of 
coordinates corresponding to tpi = ip% = 4'3 = iMo, we 
obtain from Table |TT] that the point group symmetry is 
Ds(D^) (international notation 32). That is, the chosen 
origin has a 3- fold rotation axis along [111] and three 
two-fold axes along [110], [101] and [011]. Obviously, M 
must vanish at a point of such high symmetry. Hence 
there is a point node at the origin. 

Symmetry operations consisting of a rotation followed 
by an appropriate translation yield similar point nodes 
at (with 3-fold axis along [111]), (3-fold axis 

[111]) and |i| (3-fold axis [111]). Finally, each of these 
nodes is doubled inside one unit cell because a transla- 
tion by (|,|,|) amounts to M — ► -M. The 2- and 
3-fold rotation axes form a complex array of skyrmion- 
like lines, all with winding numbers of +1. The real-space 
representation is 



M(r) = SM ( 




(23) 



; (q) 



(26) 



at momenta q = ik,^ ± k j2 ± kj 3 4£ These modes modify 
the detailed magnetic structure, but they do not change 
its symmetry, since the field respects all the symme- 
tries of the spin crystal. 

We now briefly discuss the consequences for the three 
helical magnetic structures under consideration. 

A single spin-density wave involving wavevectors ±k 
might create higher harmonics at ±3k via Eq. (|2"6")l . How- 
ever, in the case of spin spirals, = and therefore 
h,/;,3k = 0. Thus, there are no higher harmonics created 
by a single spin spiral. 

The sc spin structure with principal ordering wave- 
vectors along (001) with |k,| = Q generates higher har- 
monics along (111) (with |q| = V3Q) and along (012) 
(with |q| = VbQ). Note that throughout the current 
and last sections, all crystal directions refer to the mag- 
netic crystals. The orientation of a magnetic crystal with 
respect to the atomic crystal depends on the anisotropy 
term F a , which will be considered in Section IV A 21 

In contrast to the former cases, bec structures couple 
linearly to the q = mode (i.e., the uniform magnetiza- 
tion), since some triples of ordering vectors add to zero. 
This coupling will be further investigated in Section [VB I 
Here, we only notice that for the bccl and bcc2 states, 
h^, q =o = 0. This results can be understood in terms of 
the symmetry of these states. In the case of bccl, the 
point group symmetry is too high to support a non-zero 
Therefore, bccl and bcc2 do not 



axial vector h 



i/>,q=0- 



create a spontaneous net magnetization. The next set of 
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wavevectors which can be reached by adding three order- 
ing vectors are along (001) (with |q| = \/2Q). However 
for bccl and bcc2, direct calculation shows h^, = for 
these modes. As before, this can be understood in terms 
of symmetry. Higher harmonics along (001) would have 
the structure of a sc spin crystal. We have seen in Sec- 
tion IIV Al that the point group symmetry of sc is lower 
than that of bccl. Therefore, bccl can not create such 
an exchange field h^. We conclude that bccl creates no 
secondary Bragg peaks at (0, 0, V2Q), etc. The same is 
true for bcc2. The shortest wavevectors which are cre- 
ated by bccl or bcc2 as higher harmonics are along (112) 
(with |q| = V3Q). Others are at (110) (|q| = 2Q), (013) 
(|q| = \/5Q), (HI) (|q| = V6Q) and (123) (|q| - V7Q). 



V. RESPONSE TO CRYSTAL ANISOTROPY, 
MAGNETIC FIELD AND DISORDER. 

A. Effect of crystal anisotropy 

So far, our free energy has been completely rotation 
invariant. In the magnetically ordered states, full rota- 
tion symmetry is spontaneously broken, but any global 
rotation of the spin structure leaves the energy invariant. 
This degeneracy is lifted by an additional anisotropy term 
Fa, which couples the magnetic crystal to the underlying 
atomic lattice. The crystal anisotropy energy is small 
and may be treated as a perturbation which merely se- 
lects the directional orientation, but does not otherwise 
affect the magnetic state. 



1. Single- spiral state 

In the case of a single-spiral state, crystal anisotropy 
is a function F a (k), where k is the spiral direction. 
The function -F a (k) may depend on various parameters, 
it should be symmetric under the point group of the 
(atomic) crystal lattice and satisfy F a (k) = F a (— k). For 
concreteness, we assume the cubic point group T, rele- 
vant for the -B20 crystal structure. We further assume 
that -F a (k) is a slowly varying function, since a singular 
or rapidly oscillating function in reciprocal space would 
translate into a (non-local) interaction between magnetic 
moments and the atomic crystal. Such a function F a (k.) 
generally has its minimum at either (100) or (111), which 
can be shown in two different ways. 

The first argument is based on combining symmetry 
with Morse's theory of critical points^ Morse theory 
implies that 

maxima — saddles + minima = 2 (27) 

for a function on the unit sphere. Symmetry requires 
that F a (k) has stationary points (points with vanishing 
first derivative, i.e., maxima, minima or saddles) at (100) 



(6 directions), (111) (8 directions) and (110) (12 direc- 
tions). If F a (k) is slowly varying, we suspect that these 
are the only stationary points, since adding more max- 
ima, minima and saddles means that the function is more 
rapidly oscillating. Under this hypothesis, it follows from 
Eq. I[27p. that the (110) directions are saddle points and 
that the extrema are at (111) and (100). For (110) to be 
minima, -F a (k) needs to have additional stationary points 
(e.g. saddles) at non-symmetric, parameter-dependend 
locations. We conclude that an anisotropy which favors 
(110) would need to be more rapidly oscillating than re- 
quired by symmetry. 

The second argument is based on an expansion of F a (k) 
in powers of the directions cosines k x ,k y ,k z : 

F a (k) = a (k$ + k* + fcj) + a ' k 2 x kfk 2 z + ..., (28) 

where we retained the first two terms allowed by cubic 
symmetry. Because of the smallness of the wave-vector 
sphere Q, one typically expects \a'\ <C |a| and subse- 
quent terms even smaller. It is easily checked that for 
most values of the parameters a, a' , Eq. [55] has its global 
minima at either (100) (for a < min{0, a'/ 18}) or (100) 
(for a > max{2a'/9, a' /IS}). Only in the narrow param- 
eter regime < a < 2a'/9, the minima are indeed at 
(110)4^ We conclude that crystal anisotropy which fa- 
vors (110) may only appear in a narrow regime between 
two phases which favor (111) and (100), respectively. 

Accordingly, (111) or (100) are the selected spiral di- 
rections in all cubic helimagnets known so far.— ^ The 
preferred direction in MnSi at low pressure is (111) and 
in Fe x Coi_^Si, it is (100). In FeGe, there is a phase 
transition between these two directions^ but no inter- 
mediate phase with (110) spiral orientation has been re- 
ported. However, the neutron scattering datai^ in the 
partially ordered phase of MnSi clearly show a maxi- 
mum signal along the (110) crystal directions. While 
it is initially tempting to interpret the partially ordered 
state of MnSi as a single-spiral state that has lost it's 
orientational long-range order by some mechanism, one 
would still expect a maximal scattering intensity in the 
energetically preferred lattice direction. Theories of the 
partially ordered state in terms of disordered helical spin- 
density wave a 17 i 20 thus depend on a crystal anisotropy 
that prefers spiral directions along (110). As we have 
shown, this seems very unlikely . 

2. Helical spin crystals 

For multi-mode spin crystals, F a is no longer deter- 
mined by a single direction k, so the arguments of Section 
IV A II do not apply. Rather, the anisotropy energy de- 
pends on three Euler angles, which rotate the full three- 
dimensional magnetic structure relatively to the atomic 
crystal. In other words, F a is a function of the rotation 
group 5*0(3). Relative to some standard orientation kj 
of the mode directions, the leading-order anisotropy term 
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is 

F a (R) = a^g(Rk j )\i> j \ 2 (29) 

3 

where R is a rotation operator and gfa) = k^+ky+k^. As 
before, we have assumed a cubic point group symmetry. 

In the case a > 0, the modes of the bcc spin crystals 
get locked to the (110) directions. The orientation of 
sc is four times degenerate if a > 0. The four minima 
of F a are obtained from the standard orientation along 
(100) through a 7r/3-rotation around any of the four space 
diagonals, such that the three spiral modes point along 
(122). 

In the opposite case (a < 0), sc is oriented along (100). 
This time, it is the bcc spin crystals that get rotated by 
7r/3 around any (111) axis to reach one of four stable 
orientations. Under such 7r/3-rotations, three of the six 
modes remain along (110) and three move to (114). Each 
individual (114) direction appears only in one of the four 
solutions but each (110) direction appears in two of four 
solutions. 

If there is more than one degenerate orientation, the 
sample typically breaks up into domains such that full cu- 
bic symmetry is restored in the neutron scattering signal. 
Table UTTl lists the directions of magnetic Bragg peaks for 
the different cases. Out of the three prominent phases in 
our phase diagram, the bccl spin crystal is the only one 
that can explain the neuron-scattering peaks along (110) 
in the 'partial order' phase of MnSi. It does so most nat- 
urally for the case a > 0, which is the known sign of the 
anisotropy in MnSi at low pressure. 

In order to compare the energy scale of F a (i.e., the 
locking energy) for the different magnetic states, we note 
the following. At the phase boundary between two equal- 
amplitude spin-crystal phases (one phase with ampli- 
tudes \tp± | = ... = \iptf\ and the second phase with 
\ipi\ = . . . = \il>fj\) the amplitudes of the two neighboring 
phases are related by 

N\^ = N\^,\ 2 . (30) 

In the vicinity of the phase boundary, the anisotropy 
term is therefore proportional to the mode-average 
i/Nj^j .9(^kj). Using this result, we find that the effec- 
tive anisotropy energy is smaller for the bcc spin crystals 
than for the single-spiral state by a factor of 4 - 4.542 For 
the sc state, the locking energy is anisotropic. Certain ro- 
tations are equally costly in energy as for the single-spiral 
case while some small rotations about the minimum for 
a > are softer than for the single-spiral by a factor of 
4.5. 



B. Effect of magnetic field 

A uniform external magnetic field H couples to the 
q = mode of the magnetization, m = (M) , via Zeeman 





spiral No. 


bcc 


No. 


sc 


No. 


o > 


(Ill) 4 


(110) 


1 


(122) 


4 


a < 


(100) 3 


(110), (114)° 


4 


(100) 


1 



"When averaged over domains, Bragg peaks along (110) are twice 
as intense as peaks along (114). 



TABLE III: Crystal directions of magnetic Bragg peaks and 
number of degenerate orientations for three magnetic struc- 
tures, sc, bcc and single-spiral, with crystal anisotropy given 
by Eq. (T29J. 

coupling. The uniform magnetization, in turn, couples 
to the helical modes tpj through 

F = F \ m =o + r o m 2 + 17(0, 0, 0) m 4 - h^O) m 

+ 2 (Ufa, -k is 0) m 2 + Ufa, 0, 0) m 2 ± J) |^| 2 , 

(31) 

where mj_j = m (m • kj)kj and we have used 
Eqs. ([5^1.(p5]). Plumer and Walker^ 8 , argued that 
17(0,0,0) « U(kj,-kj,0) « Us, which we will use in 
the following for simplicity. 

1. Response of the single-spiral state 

The behavior of the single-mode spiral state 
under a magnetic field has been studied both 
exp er iment ally&Sai i 1 1 1 2 and theoretically 37 ' 38 It is char- 
acterized by a strongly anisotropic susceptibility, induced 
by the last term in Eq. (|31[) . For a fixed spiral direction k, 
the susceptibilities parallel and orthogonal to the spiral 
direction are given by 

/ U' V 1 (32) 
x± «(A + 2— |r(Q)|J , 

where A = 2[r - r{Q)\ = 2JQ 2 and U' = J7(k,0,0). 
Well below the critical ordering temperature, A <C 
2|r(Q)| and therefore \\\ ^ X±- This strong anisotropy 
leads to a spin reorientation transition at H — H sr , where 
the spiral axis gets oriented along the field direction. The 
value of H SI depends on the anisotropy [Eq. ([28]) ] and the 
field direction. For a > and H || (100), Plumer and 
Walker obtained 

4rv 

Hi = — > 4aA, (33) 

X|| - X-L 

where we have used x\\ 3> X-L > 0- Once the spiral is 
oriented, the susceptibility is large (equal to x\\)- 

The spiral amplitude decreases as a function of the 
external field and vanishes at 

H c = |Vo|A, (34) 

where IV'ol 2 = \ r {Q)\/(2U s ). Above H c , the magnetiza- 
tion is uniform. 
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2. Response of the bccl spin crystal 

In the bec spin crystal states, the linear response is 
isotropic, because their symmetry group does not al- 
low for an anisotropic susceptibility tensor. As a con- 
sequence, there is no orientation of the bec state towards 
the magnetic field at the level of linear response (i.e., from 
energies up to order H 2 ). However, there is a sub-leading 
contribution to the energy oc {M x M y M z )H x H y H z . This 
contribution splits the degeneracy between the S = 1 and 
S = — 1 states and it may lead to a reorientation of the 
bec crystal towards the field. 

In terms of the six ^-variables of bec, the exchange 
field h^(0), which enters Eq. (|31[) . amounts to 



M°) = -juRe[(5-n/2)h 



(35) 



where fi = U(kj 1 , kj 2 , kf 3 ) for kj 1 + kj 2 + kj 3 = and 



mia \ i J - ami ( -j 



(36) 



The (isotropic) inverse spin susceptibility (see Ap- 
pendix) in the bccl state is composed of three contri- 
butions 



-^bccl 

The first term 



•^barc 



Y _1 
phase 



Xbarc = 2 r 



ljj> 
3 I 



<Q)\ 



(37) 



(38) 



where U hccl = 1/6[U S + V p (tt/4) + 4V p (w/6) - A nt ], can 
be derived in analogy to the single-spiral case. In fact, 
Xbarc is a "mixture" of xil an d X-Lj determined geometri- 
cally by the angles between the mode directions k, and 
the magnetic field. It follows that Xbare "C Xll , provided 
C^bcci ~ U s (the two couplings are equal at the phase 
boundary between single-spiral and bccl). The remain- 
ing terms in Eq. (f3"T)) 



m 2 KQ)I 



3 C^bcci An 



X. 



25m 2 |KQ)I 



(39) 



amp 



12£/ b cci [U s -V p (ir/4) + \ n 



stem from the response of the bec magnetic structure to 
the field. That is, they originate from the adjustments of 
relative phases and amplitudes, respectively, of the heli- 
cal modes as a result of the term — h^,(0) • m in Eq. (f3"Tj) . 
The effect of 

Xphase and Xampj which are necessarily neg- 
ative, is to increase the susceptibility of the bccl state. 

The change in the relative amplitudes and phases of 
the six interfering spirals as a function of the magnetic 



field may be calculated (see Appendix). For example, the 
linear response of the amplitudes of bccl is 



0~\i>2 

S\A 

V^ 6 



\ 



5 fiS 



4p7 a - V p (tt/A) + A nt ] 



— m 
rn 

—m 
rn 
m 



\ 



X 



X 



(40) 



V -my J 



This response should be observable by neutron scatter- 
ing, if it is possible to prepare the sample in a single- 
domain state (i.e. without mixture of the two time- 
reversal partners). For example, a field in z direction 
affects \ipi\ and l^l, the amplitudes of the modes prop- 
agating orthogonally to z (Fig. [5]), which get enhanced 
and suppressed by the magnetic field, respectively. 

The expected effects of external magnetic field on the 
resistivity of bec spin crystal are presented elsewhere^ 



VI. EFFECT OF IMPURITIES: A POSSIBLE 
ROUTE TO 'PARTIAL ORDER'. 

While the helical spin crystal states are expected to 
show Bragg spots at particular wavevectors, a variety of 
effects such as thermal or quantum fluctuations or dis- 
order can destroy the long range order while preserving 
the helical spin crystal structure at shorter scales. Here, 
we investigate in more detail the effect of non-magnetic 
disorder on helical spin crystal structures. 

Although the experimentally studied helimagnets are 
very clean from the electrical resistivity point of view, the 
helical magnetic structures are sensitive to disorder at a 
much longer length scale. In addition, the low energy 
scales required to distort them means that one needs to 
consider disorder effects. An observation that can imme- 
diately be made is that for the physically relevant case 
of non-magnetic disorder (V,u s (r)), the single-spiral state 
and the spin crystal states respond very differently By 
symmetry, the coupling of disorder to the magnetic struc- 
ture is given by F dls = (Vd is (r)|M(r)| 2 ). Hence, single- 
spiral states which are unique in having a spatially uni- 
form magnitude of magnetization (|M(r)| = constant) 
are unaffected by this coupling; in contrast the spin- 
crystal states necessarily have a modulated magnitude 23 
and hence are affected by non-magnetic disorder. There- 
fore the neutron scattering signal of the spin-crystal state 
is expected to have more diffuse scattering than the sin- 
gle mode state. This is consistent with the experimental 
observation that the high pressure phase has diffuse scat- 
tering peaked about (110) while the low pressure phase 
has sharper spots, consistent with identifying the two as 
spin-crystal and single-spiral states respectively. 

The effect of disorder on the spin-crystal state is closely 
related to the problem of the ordering of an XY model 
in the presence of a random external field. The phase 
rotation symmetry of the XY model captures the trans- 
lational invariance of the spin-crystal in the clean state. 
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Disorder destroys this invariance and behaves like a ran- 
dom field applied to the XY system. Using the insights 
from the study of that problem in three dimensions,— one 
expects that for weak disorder a Bragg glass will result, 
where although true long range order is destroyed, power 
law divergent peaks at the Bragg wavevectors remain, 
and the elastic constants remain finite. For stronger 
disorder one expects this algebraic phase to also be de- 
stroyed, and recover a short range correlated phase with- 
out elasticity. Nevertheless, for the case of the beel and 
bcc2 crystals, due to time reversal symmetry (T) break- 
ing in these states, the disordered states also sponta- 
neously break time reversal symmetry, and hence a phase 
transition is expected on cooling despite the absence of 
long range order. It is difficult to predict which of these 
two scenarios (Bragg glass or only T breaking) is more 
appropriate for MnSi. In the latter case one may esti- 
mate the spreading of the Bragg spots due to disorder by 
considering the energetic cost to deform the spin-crystal 
state in different ways. 

Ignoring elastic contributions, there are two distinct 
types of deformations - ones that involve a change in the 
magnitude of the ordering wavevectors 5qu and others 
that do not change the wavevector magnitude but rotate 
the structure from its preferred orientation: <5qj_. The 
second is expected to be low in energy because rotations 
of the structure are locked by the crystal anisotropy term, 
which is weak. From Eq. (|28|) . we obtain the energy 
cost to shift the ordering vector by Sq± along the sphere 
|q| = Q 



8F ± 



4a / dq± 



3k \ Q 



(41) 



where k = 1 for the single-spiral. For multi-mode spin 
crystals, the energy cost of rotation is reduced, as ex- 
plained in Section IV A 21 Thus, k 4 for the bec spin 
crystals. In contrast, deformations that change the mag- 
nitude of the ordering wavevectors, must contend with 
the DM interaction scale, and hence pay a higher energy 
penalty 



<5F|| = -A -l^o | 



Sq\\ 
Q 



(42) 



Assuming the disorder couples to these deformations 
equally, we can estimate the ratio of their amplitudes 
in the limit of weak deformations, by equating Eqs. (fITj) 
and ([32). It follows 



Sq± 
&?|| 



3kA|V> | 
8a 



(43) 



Using Eqs. (j3"3"|) and ((3~4")) , we can relate this ratio to the 
experimentally known ratio between the critical mag- 
netic fields for, respectively, reorienting and polarizing 
the single-spiral state 



5q± > 3k H c 



(44) 



We can now apply these results to the case of MnSi 
and test the hypothesis that the 'partial order' state 
is in fact a disordered bec spin crystal. Setting in 
k = 4 and the experimentally measured^ critical fields 
for MnSi, H c = 0.6Tesla and H c i — O.lTesla, one obtains 
5q±/5q\\ > 15. Neutron scattering experiments do indeed 
find that the transverse broadening is larger than the lon- 
gitudinal broadening, but since the latter is resolution 
limited, this only gives us an lower bound that is consis- 
tent with the estimate above: [<5g_i_/<5<?||]cxpt > 2.3. Nev- 
ertheless, the trend that the width of the spot is greater 
along the equal magnitude sphere than transverse to it 
is clearly seen in the experimental data. 

Thus, weak non-magnetic disorder of the atomic crys- 
tal is expected to destroy magnetic long range order in 
multi-mode helical spin crystal states and lead to a neu- 
tron scattering signal compatible with the observations 
in the 'partial order' phase of MnSi. However in the case 
of bec spin crystals, time-reversal symmetry breaking is 
expected to persist even in the presence of disorder. The 
scenario of interpreting 'partial order' in MnSi as a beel 
state disordered by impurities thus predicts quasi-static 
local magnetic moments and implies a finite temperature 
phase transition on cooling into this phase. 



VII. CONCLUSION 

We have analyzed the magnetic properties of non- 
ccntrosymmctric weak ferromagnets subject to DM spin- 
orbit coupling. This problem falls into the general class 
of systems where the low energy excitations live on a 
surface in reciprocal space rather than on discrete points. 
The addition of DM interactions to a ferromagnetic state 
produces a large degeneracy of magnetic states char- 
acterized by arbitrary superpositions of spin helices of 
a fixed helicity and fixed wavevector magnitude. This 
enormous degeneracy is broken by interactions between 
modes, and the single-spiral state is realized for slowly 
varying interactions, by virtue of its unique property of 
having a spatially uniform magnitude of magnetization. 
For more general interactions, multi-mode helical spin 
crystal states are obtained. We show that for the model 
interactions considered, the phase diagram is largely de- 
termined just by considering the interactions between 
pairs of modes. The phase that is eventually realized 
may be readily deduced from the range of angles in which 
this interaction drops below a critical value. In particu- 
lar, the bec structure is stabilized by virtue of the fact 
that its reciprocal lattice, fee, is a close packed structure. 
These results may also be relevant in other physical sit- 
uations where crystallization occurs, such as the Larkin- 
Ovchinnikov-Fulde-Ferrel instability in spin-imbalanced 
superconductors, which may potentially be realized in 
solid state systems^ cold atomic gases^i and dense nu- 
clear matter 42 

Helical spin crystals typically give rise to complicated 
real space magnetic structures which we discussed in this 
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paper. In particular, topological textures like merons and 
anti-vortices can be seen about special axes in particu- 
lar realizations, although these are not expected to be 
stable given the order parameter and spatial dimension- 
ality of the system. We show here that such topological 
structures exist as a consequence of symmetry, which also 
dictates the absence of certain higher Bragg reflections, 
which a naive analysis would predict. 

The response of helical spin crystals to crystalline 
anisotropy and applied magnetic field are considered with 
a special emphasis on the bcc structures which are con- 
trasted against the response of the single helix state. An 
unusual transfer of spectral intensity in the presence of 
an applied magnetic field, which is strongly dependent 
on the direction of applied field is noted for the bcc 
structures. This is a consequence of broken time reversal 
symmetry in the absence of a net magnetization (which 
is symmetry forbidden). The unusual magnetotransport 
in such a state, a linear in field magnetorcsistance and 
quadratic Hall effect, has been discussed briefly in Rcf. 23 
and was elaborated upon in Refill. 

Helical spin crystals exhibit Bragg peaks at specific 
wave- vectors, and hence are not directly consistent with 
the experimental observation of 'partial order'. The point 
of view taken in our earlier work2£ is that the short 
distance and short time properties are captured by the 
appropriate helical spin crystal structure. Studying the 
properties of helical spin crystals with long range order is 
a theoretically well defined task with direct consequences 
for a proximate disordered phase with similar correlations 
up to some intermediate scale. The mechanism that leads 
to the destruction of long range helical spin crystal order 
is uncear; in Refi^, this was assumed to be the coupling 
to non-magnetic disorder. Then, as elaborated in this pa- 
per, beginning with a bcc helical spin crystal a neutron 
scattering signature consistent with that of 'partial-order' 
may be obtained. However, within the simplest version 
of this scenario, one also expects a finite temperature 
phase transition where time reversal symmetry breaking 
develops, and static magnetic order which may be seen 
in nuclear magnetic resonance or muon spin rotation ex- 
periments. Other mechanism for the destruction of long 
range order of the bcc spin crystal state, such as thermal 
or quantum fluctuations may also be considered, but are 
left for future work. 



Acknowledgments 

We would like to thank L. Balents, P. Fazekas, I. Fis- 
cher, D. Huse, J. Moore, N. Nagaosa, M.P. Ong, C. Pflci- 
derer, D. Podolsky, A. Rosch, T. Senthil, H. Tsunetsugu 
and C. Varma for useful and stimulating discussions and 
V. Aji for an earlier collaboration on related topics. This 
work was supported by the Swiss National Science Foun- 
dation, the A. P. Sloan Foundation and Grant No. DE- 
AC02-05CH11231. 

APPENDIX: LINEAR RESPONSE 

Let the free energy F depend on real internal variables 
x = (x\, . . . ,x n ) (related to various order parameters) 
and on m = (mi, 7712,777,3), which couples linearly to the 
external field H. The internal variables are chosen such 
that for m = 0, the minimum energy shall be at x = 0. 
Expanding F to second order in x and m yields 

F = ^x T Ax + m T Bx + im T Cm, (A.l) 

where A, B and C are matrices, A = A T , C = C T . All 
eigenvalues of A are positive. Eq. (|A. 1|) can be written 
as 

F = i(x - x m ) T A (x - x m ) + imV'm. (A.2) 

where x m = -A~ 1 B T m and x" 1 = C - BA~ 1 B T . It 
follows that under the influence of an external field H, 
the equilibrium internal variables get shifted to x = x m 
and the linear response is given by m = ^H. Thus, 
there are two contributions to the inverse suszeptibility: 
Xbarc = C an d Xx 1 = ~BA~ l B T . The latter comes from 
the internal response of the x-variables to the magnetic 
field. 

These general results are applied in Section IV B 21 
where the internal variables x%, . . . ,xq are the deviations 
from their equilibrium value at m = of six amplitudes 
and three phases (holding the other three phases 
fixed). In this case, \x leads to both the phase and am- 
plitude related terms in Eq. Q39p, 
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Whenever W' < 0, Eq. 1)160 becomes negative for a rapidly 
varying M 2 and the quartic theory is thus unstable. We 
consciously ignore this instability, since we only study local 
minima at small M 2 , where F? is dominant and forbids 
rapid variations. 

Only the relative orientation between these modes enters 
the calculation, not their orientation with respect to the 
atomic lattice, since crystal anisotropy is neglected. 
To obtain the phase diagram of Figs. [2] and [3] we have used 
the approximation (l + z) 5/2 + (1 -xf /2 » 2 + (4 v / 2- 2)x 2 
for x = cos 6 in the evaluation of V p (8). The error is smaller 
than 0.7%. 

The shift to directions [±y/l + S ±s/l + Sy/1- 25] is very 
small with S ~ 0.02 and even without the deformation, fee 
would still beat the neighboring phases in the fee* -region. 
In non-magnetic crystals, the coupling term is cubic and 
therefore higher harmonics are generated at all sums of 
two momenta ± k J2 . In magnetic structures, higher- 
harmonics wavevectors are restricted to sums of three or- 
dering momenta. 

In this regime, all 26 high-symmetry points are ex- 
trema of F a and 24 saddle points are located at 
(VZctVZaVa' — 4a), in agreement with Eq. (f27|) . 
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Using max{F a } — min{F a } to estimate the locking en- 
ergy scale leads to a reduction factor of 4.5. Expanding 
Fa around its minimum for a > leads to a factor of 4. 



